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Foreword 



The present document describes the detailed mapping of the general audio service employing the aacPlus general audio 
codec within the 3GPP system. 

The contents of the present document are subject to continuing work within the TSG and may change following formal 
TSG approval. Should the TSG modify the contents of this TS, it will be re-released by the TSG with an identifying 
change of release date and an increase in version number as follows: 

Version x.y.z 

where: 

X the first digit: 

1 presented to TSG for information; 

2 presented to TSG for approval; 

3 Indicates TSG approved document under change control. 

y the second digit is incremented for all changes of substance, i.e. technical enhancements, corrections, 
updates, etc. 

z the third digit is incremented when editorial only changes have been incorporated in the specification; 
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Scope 



This Telecommunication Standard (TS) describes the Parametric Stereo encoder part of the Enhanced aacPlus general 
audio codec [4] . 



2 Normative references 

This TS incorporates by dated and undated reference, provisions from other publications. These normative references 
are cited in the appropriate places in the text and the publications are listed hereafter. For dated references, subsequent 
amendments to or revisions of any of these publications apply to this TS only when incorporated in it by amendment or 
revision. For undated references, the latest edition of the publication referred to applies. 

[1] ISO/IEC 14496-3:2001/AMD1:2003: "Bandwidth Extension". 

[2] ISO/IEC 14496-3:2001/Amd.l:2003/DCOR1. 

[3] ISO/IEC 14496-3:2001/ Amd.2:2004: "Parametric Coding for High Quality Audio". 

[4] 3GPP TS 26.401: "Enhanced aacPlus general audio codec; General Description:. 

3 Definitions, symbols and abbreviations 

3.1 Definitions 

For the purposes of this TS, the following definitions apply: 

hybrid QMF: a QMF filterbank combined with additional filters to achieve higher frequency resolution for the 
lower QMF bands 

stereo band: a group of consecutive hybrid QMF subbands used for coding one stereo parameter 

3.2 Symbols 

For the purposes of this TS, the following symbols apply: 

/^ («) Subsample in hybrid QMF matrix: left channel, band k, subsample n. 

r^^ («) Subsample in hybrid QMF matrix: right channel, band k, subsample n. 

3.3 Abbreviations 

For the purposes of this TS, the following abbreviations apply. 

SBR Spectral Band Replication 

AAC Advanced Audio Coding 

aacPlus Combination of MPEG-4 AAC and MPEG-4 Bandwidth extension (SBR) 

Enhanced aacPlus Combination of MPEG-4 AAC, MPEG-4 Bandwidth extension (SBR) and MPEG-4 

Parametric Stereo 
QMF Quadrature Mirror Filter 

MPEG Moving Picture Expert Group 

IID Inter Intensity Difference, (stereo parameter) 

ICC Inter Channel Coherence, (stereo parameter) 
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Outline description 



This TS is structured as follows: 

Section 5.2 describes the hybrid QMF filterbank and its integration in the Parametric Stereo system. 

Section 5.3 describes the hybrid QMF filterbank and its integration in the Parametric Stereo system. 

Section 5.4 describes the parameter estimation algorithms and quantization. 

Section 5.5 describes how to convey the estimated parameters in the bitstream. 

Section 5.6 and section 5.7 describes preparation of the signal that should feed the aacPlus mono encoder after the 
Parametric Stereo encoding. 



Parametric stereo encoder 



5.1 System overview 





L 
























AAC 
bitstream 


X 

Z) 

E 
a 

w 
5 


■^~° 
5 « 

03 ^ 
2 CO 

0-' 




0) 


L ^ 


o 

c 
o 
E 
o 

X 

E 
c 
S 

o 
Q 


M 


II 

■a CO 
So 


M 


0) 
CD w 

2 2. 
O 


M 

► 




"D 
O 
O 

c 



o 
< 
< 


R 




R 










SBR 
bitstream 










































CD 

■a 

8 

c 

CD 
CC 

m 














PS 
bitstream 






Stereo 
parameter 
extraction 


iid(b), 
icc(b) 
















PS 
bitstream 
formatting 






^ 





































Figure 1 : Encoder overview 

The interface between the parametric stereo encoder tool and the aacPlus encoder is depicted in Figure 1 . In the figure L 
and R denotes the left and right channel respectively, while M denotes the down-mixed mono signal which the aacPlus 
encoder operates on. 

The parametric stereo coding tool is able to capture the stereo image into a limited number of parameters, requiring only 
a small overhead of a few kbit/s. Together with a controlled monaural downmix of the stereo input signal, the 
parametric stereo coding tool is able to regenerate the stereo signal at the decoder side. 

The encoder operates as a non-modifying analyzer prior to the aacPlus encoder, though it shares the same QMF analysis 
filterbank. The decoder operates as a post process to aacPlus using the Parametric Stereo data conveyed by the bitstream 
to synthesize the stereo properties of the output signal. Part from the parametric stereo tool, the aacPlus runs in mono 
mode not affected by Parametric Stereo. 

The bitstream syntax and decoder description of the parametric stereo tool in combination with aacPlus is defined in [3]. 
This system includes only the baseline level defined in that standard. 
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5.2 Analysis filterbank 

5.2.1 QMF analysis filterbank 

This filterbank is identical to the 64 complex QMF analysis filterbank as defined in ISO/IEC 14496-3/AMDl:2003, sub 
clause 4.B.18.2 [1], [2]. However, in the equation for matrix M(k,n) and in Figure 4.B.20, the term "(2*n+l)" has to be 
substituted by "(2*n-l)". The input to the filterbank are blocks of 64 samples of the monaural synthesized signal M. For 
each block the filterbank outputs one slot of 64 QMF samples. 

5.2.2 Low frequency filtering 

The lower QMF subbands are further split in order to obtain a higher frequency resolution enabling a proper stereo 
analysis and synthesis for the lower frequencies. To achieve those, in total 77 frequency bands, a hybrid filterbank 

configurations have been defined. The filter used for this sub subband filtering, Q^ is defined according to: 



G^=g'{n)exp 



.In 



r + i|(«-6) 



where g '' represents the prototype filters in QMF subband p. Q'' represents the number of sub-subbands in QMF 
subband p, q the sub-subband index in QMF channel p and n the time index. The prototype filters are all of length 13 
and have a delay of 6 QMF samples. The prototype filters are listed in Table 1. 

Table 1 : Prototype filter coefficients for the filters that split the lower QMF subbands 



n 


g''{n),Q'=i 


g'Hn),Q^'-=4 





0.00746082949812 


-0.00305151927305 


1 


0.02270420949825 


-0.00794862316203 


2 


0.04546865930473 





3 


0.07266113929591 


0.04318924038756 


4 


0.09885108575264 


0.12542448210445 


5 


0.11793710567217 


0.21227807049160 


6 


0.125 


0.25000000000000 


7 


0.11793710567217 


0.21227807049160 


8 


0.09885108575264 


0.12542448210445 


9 


0.07266113929591 


0.04318924038756 


10 


0.04546865930473 





11 


0.02270420949825 


-0.00794862316203 


12 


0.00746082949812 


-0.00305151927305 



Figure 2 and Figure 3 illustrate the hybrid analysis and synthesis filterbank for the 77 frequency bands configuration. 
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Figure 2: Hybrid QIVIF analysis filterbank providing 77 output bands. The three lower subbands of the 
64 QMF (see dashed box) are further split to provide for increased resolution for the lower 

frequencies 
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Figure 3: Hybrid QIUIF synthesis filterbanlt using 77 input bands. The coefficients offering higher 
resolution for the lower QMF subbands are simply added prior to the synthesis with the 64 subbands 

QMF (see dashed box) 

In order to time align all the samples originating from the hybrid filterbank, the remaining QMF subbands that have not 

been filtered are delay compensated. This delay amounts to 6 QMF subband samples. This means Gg (z) = Z~ for 

k=3...63. In order to compensate for the overall delay of the hybrid analysis filterbank, the first 10 sets (6 from delay 
and 4 from QMF filter) of hybrid subbands are flushed and therefore not taken into account for processing. 

The resultant of this operation is a slot of hybrid subband samples consisting of a LF (low frequency) sub QMF subband 
portion and HF (high frequency) QMF subband portion. 



5.3 Configurations 



The parametric stereo encoder uses two different configurations depending on desired frequency resolution. The 
configuration parameter, num_stereo_bands determines what frequency resolution should be used for the stereo 
parameters. For all bitrates below 21000 bit/s, num_stereo_bands is set to 10 otherwise num_stereo_bands is set to 20. 



5.4 Stereo parameter extraction 



5.4.1 Parameter estimation 

In order to estimate the stereo parameters the signals M, L and R are analyzed using the hybrid filterbank as in Figure 2 
for providing the 77 frequency bands addressed by the index, 0<k <11 . This results in the (sub-)subband domain 

signals, m[k,n) , l[k,n) and r{k,n) . 

To estimate the parameters for the current frame the following is calculated: 



h L-i 

k=ki M=0 



l\ k, + l + « 

2 



+ £ 
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** L-l 



-.(^)=ZZ 



r\ k, + l + « 

2 



+ e 






/ 



« r 



\ 



A:, + l + n 

V 2 



+ £• 



7 



where e, (Z?) , e^ (Z?) and e^ (Z?) are the left channel excitation, the right channel excitation and the non-normalized 

cross-channel excitation between left and right channel for stereo bin b respectively, L the segment length, S a very 

small value preventing division by zero {£ = 10" ). The summation over k is shown in Table 2 for 20 bands case. For 
the 10 bands case additional formulas below are used. 

Table 2: Summation range in 71 sub subbands in case of 20 bands 



Parameter index b 


Subsubband index 


QMF cinannel 











1 


1 





2 


2 





3 


3 





4 


10 


1 


5 


11 


1 


6 


12 


2 


7 


13 


2 


8 


16 


3 


9 


17 


4 


10 


18 


5 


11 


19 


6 


12 


20 


7 


13 


21 


8 


14 


22-23 


9-10 


15 


24-26 


11-13 


16 


27-30 


14-17 


17 


31-35 


18-22 


18 


36-47 


23-34 


19 


48-76 


35-63 



For the 10 band case, summation will be use the same table as the 20 band case but with the additional summation: 

e^{2b) + e^{2b + \) 



%o(^) = - 



2^,10 (^) = 



««,io(^) = 



e^{2b) + e^{2b + l) 



e,{2b) + e,{2b + l) 



Where e, j^ and e^,Q and e^,Q should replace e, and e^ and e^ in the following formulas for the 10 band case. 
The IID, denoted as iid(b) and the ICC, denoted as icc(b), for each stereo band b are calculated as: 
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iid {b) = lOlog^ 






ICC 



{b) = 






1-min 



^e,{b)e^{b) 



f 



1-min 



'Ab) 



h^b)eXb) 



,1 



,h<5 and num _ stereo _bands =10 or 
,b<ll and num _ stereo _bands = 20 

, otherwise 



5.4.2 Quantization of IID and ICC parameters 

The obtained values for IID and ICC are quantized to the nearest lower value given in Table 3 and Table 4, respectively, 
then coded into Huffman words according to Huffman tables found in [4] section 8.B.1 Huffman tables ps_data. 

Table 3: Quantization grid for iid 



Index 


-7 


-6 


-5 


-4 


-3 


-2 


-1 





IID [dB] 


-21.5 


-16 


-12 


-8.5 


-5.5 


-3 


-1 





Index 


1 


2 


3 


4 


5 


6 


7 




IID [dB] 


1 


3 


5.5 


8.5 


12 


16 


21.5 





Table 1 : Quantization grid for ice 



index 





1 


2 


3 


4 


5 


6 


7 


Ice 





0.088900 


0.229800 


0.364250 


0.504500 


0.635100 


0.799600f 


0.945650 



5.5 Writing to bitstream 



The parametric stereo bitstream should be placed in the SBR extension field according to syntax in Table 5. The 
EXTENS10N_1D_PS should be used as extension identifier. Note that due to SBR extension is design, the total size, in 
bytes, of the ps_data element has to be added to the SBR extension element before writing the actual ps_data. 

Table 5: Syntax of sbr_extension() 



Syntax 



No. of bits Mnemonic 



sbr_extension(bs_extension_id, num_bits_left) 

{ 

switcin (bs_extension_id) { 
case EXTENSION_ID_PS: 

num_bits_left -= ps_data(); 

break; 
default: 

bs_fill_bits; 

num_bits_left = 0; 

break; 



Note 1 



num bits left bslbf 



± 



} 



Note 1 : ps_data() returns the number of bits read. 
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Table 6: Values of the bs extension id field 



Symbol 


Value 


Purpose 


EXTENSION_ID_PS 


2 


Parametric Stereo Coding 




all other values 


reserved 



The syntax of the ps_data is defined in [4] Table 8.1. The bitstream elements in ps_data should be assigned values 
according to the following list. 

enable_ps_header - Is set to 1 for all frames containing an SBR header. Is set to for all other cases. 

enable_iid - Is set to if all values in the iid vector after quantization are zero. Is set to 1 for all other cases. 

enable_icc - Is set to if all values in the ice vector after quantization are zero. Is set to 1 for all other cases. 

iid_mode - Is set to if 10 stereo bands resolution is chosen. Is set to 1 for all other cases. 

icc_mode - Is set to if 10 stereo bands resolution is chosen. Is set to 1 for all other cases. 

enable_ext - Is set to 0. 

frame_class - Is set to 0. 

num_env_idx - Is set to if all values in the iid and ice vectors after quantization are unchanged since previous frame. 
Is set to 1 for all other cases. 

iid_dt - Is set to 1 if differential coding of iid parameters over time gives a lower total Huffman code length than 
differential coding of iid parameters with respect to the previous parameter position. Is set to for all other cases. 

icc_dt - Is set to 1 if differential coding of ice parameters over time gives a lower total Huffman code length than 
differential coding of iid parameters with respect to the previous parameter position. Is set to for all other cases. 

iid_par_dt[n] - In case of differential coding of IID parameters over time (iid_dt==l), iid_par_dt[n] describes the IID 
difference with respect to the previous parameter position. If no previous parameter position is available, iid_par_dt[n] 
represents the IID difference with respect to the decoded value (i.e. index=0). 

iid_par_df[n] - In case of differential coding of IID parameters over frequency (iid_dt==0), iid_par_df [n] describes the 
Huffman encoded IID difference with respect to the parameter (n-1). 

iid_par_df[0] - represents the IID difference with respect to the decoded value (i.e. index=0). 

icc_par_dt[n] - In case of differential coding of ICC parameters over time (icc_dt==l), icc_par_dt[n] describes the 
difference with respect to the previous parameter position. If no previous parameter position is available, icc_par_dt[n] 
represents the ICC difference with respect to the decoded value 1 (i.e. index=0). 

icc_par_df[n] - In case of differential coding of ICC parameters over frequency (icc_dt==0), icc_par_df[n] describes 
the ICC difference with respect to the parameter (n-1). icc_par_df[0] represents the ICC difference with respect to the 
decoded value 1 (i.e. index=0). 



5.6 Downmixing to mono 



A monaural signal is needed for encoding of the HE AAC mono core. A simple stereo-to-mono downmixing scheme is 
performed after extraction of the stereo parameters. It combines equal parts of left and right channel and scales the 
resulting mono signal in order to preserve most of the original total power. The process is defined according to: 

, , l(k,n) + r{k,n) , , 
where the stereo scale factor, y is defined by: 
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Y[k,n) = min 



2,. 



l[k,n) + r[k,n) 



'0.5-\l{k,n) + r{k,n)\' 



0<k<NUM OF BANDS 



5.7 Synthesis filterbank 



The stereo to mono down-mixed hybrid subband signal m{k,n) is fed into the hybrid synthesis filterbank, that is 
implemented as an adder of sub QMF samples. This is illustrated in Figure 3. 

The synthesis filtering and implicit down-sampling of the 64 subband signals is achieved using a 32-channel QMF 
bank. The output from the filterbank is real-valued time domain samples. The process is given by the flowchart in 
Figure 4. The synthesis filtering comprises the following steps, where an array v consisting of 640 samples is assumed: 

• Shift the samples in the array v by 64 positions. The oldest 64 samples are discarded. 

• The array of 32 complex-valued subband samples Z is separated into the real and imaginary components as Z = R 
+ i I. The components are scaled and DCT and DST type IV transformed as 



1 '' 
r{n) = —Y,R(h) COS. 

1 ^' 



—{k + -){n + -) 

_32 2 2 

— {k + -){n + -) 

32 2 2 



,0<«<32 



• The arrays r and i are combined and stored in the positions to 63 in array v as 

v{n) = i{n)-r{n) ,0<?i<32 

v{6?>-n) = i{n) + r{n) ,0<?i<32 

• Extract samples from V according to the flowchart in Figure 4 to create the 320-element array g. 

• Multiply the samples of array g by every other coefficient of window w. The window coefficients of C can be 
found in Figure 4, and are the same as for the analysis filterbank. 

• Calculate 32 new output samples by summation of samples from array w according to the last step in the flowchart 
of Figure 4 

Every SBR frame produces an output of numTimeSlots ■ RATE ■ 32 time domain samples. In the flowchart of Figure 
4 X[k][l] corresponds to subband sample 1 in the QMF subband k, and every new loop produces 32 time domain 
samples as output. 
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Start 
[ for QMF subsample I 



for( n = 639; n >= 64; n--) { 

v[n] = v[n -64] 
} 



for( n = 0; n< 32; n++) { 
31 

r[n] = 1/ 64 * Z R[k][l] * cos( :t / 32 * (k + 0.5) * (n + 0.5) 

k=0 
31 

i[n] = 1/ 64 * Z l[k][l] * sin( n / 32 * (k + 0.5) * (n + 0.5) 
k=0 



for 


(n = 0; 


n < 


32; 


n++){ 




v[n] 


= 


i[n] 


-r[n] 




v[63-n 


= 


i[n] 


+ r[n] 


} 











for( n = 0; n <= 4; n++) { 
for( k = 0; k <= 31 ; k++) { 

g[64*n + k] = v[128*n + k] 

g[64*n + 32 + k] = v[1 28 * n + 96 + k] 



} 



} 



for( n = 0; n <=319; n-i-i-) { 

w[n] = g[n] * c[2*n] 
} 


1 


r 



for( k = 0; k <= 31 ; k++) { 
temp = w[k] 

for( n = 1 ; n <= 9; n++) { 
temp = temp + w[32*n + k] 

} 

nextOutputAudioSample = temp 

} 



Done 



Figure 4: Flowchart of encoder synthesis QMF bank 
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